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1 • Introduction 


We consider the initial value problem for a system of conservation laws 
given by 


( 1 . 1 ) 




+ 




u, f € If 1 


u(0 , x 0 ,***,x d ) = <J 

( u _( x 0 » ••*> x d )» 


x 0 > S(x 1 ,»»»,x d ) 
x 0 < S(x 1 ,***,x d ) 


and satisfying 

(1) f^, S are analytic, u + , u_ are analytic across S; however, 
u (0 , x ,— ,x ) may have a small jump discontinuity or a jump in 
derivatives, not necessarily small, at S. 

(2) Equation (1,1) is hyperbolic in the following sense: If we let 

d af . 

M(uj,u) = £ a) -r— , a) 6 - {0}, then M has real eigen- 

i«0 1 9u 

values X^(o),u) _< X^C^u) _< ,•••, _< X n (o),u) and a basis of eigen- 
vectors r^ (o),u) , • • • ,r^((o,u) . We denote left eigenvectors by £.(o),u). 

(3) Equation (1.1) has either genuinely nonlinear or linearly degenerate 
fields, i.e., 


either 


or 


V X. *r . * 0 
u 1 1 


V X «r . =0 
u i l 


for u in a neighborhood of u(0,Xq , • • • ,x^) and |u)| = 1. 



2 


(4) If M has multiple eigenvalues, then the corresponding field must be 
linearly degenerate. 

Our object is to obtain a power series representing a distribution 
solution to (1.1). 

The conditions (2), (3), (4) are in part dictated by the properties 
characterizing the Euler equations. For a polytropic three-dimensional gas 
flow they are given by 


— — 


— — 




— — 


— i 

p 


pu 


pv 


pw 


0 

pu 


2 A 

pu + p 


puv 


puw 


0 

pv 

+ 

puv 

+ 

2 

pv + p 

+ 

pvw 

= 

0 

pw 


puw 


pvw 


2 , 

pw + p 


0 

e 

t 

u(e + p) 

X 

v(e + p) 

y 

w(e + p) 

z 

0 


9 9 9 

with p = (y-l)(e-(p(u + v + w )/2)) where p = density, u,v,w 

velocities, e = total energy, and p = pressure [2]. 

The eigenvalues of M, in this case, are 


(DqU + (OjV + U^W - C < 0 )qU + U)jV + (O^W < <i)qU + (DjV + u^w + c 


where c = /4^- is the speed of sound in the medium. The first and third 
/dp 
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fields are genuinely nonlinear and the corresponding eigenvalues simple. The 
second field is linearly degenerate with the eigenvalue of multiplicity 3. 
There is, however, a basis of eigenvectors so (2), (3), (4) are satisfied. 

As a preliminary step we change variables to make the initial surface of 
discontinuity flat. If 


X = X Q - S(x lf -*«,x d ) 


y 


i 


x. 

l 


i = 1,2, 


>d 


t = t, 


then from (1.1) 


( 1 . 2 ) 


u t + ( f 0 (u ’ y) )x + X 
1=1 


u(0 ,x,y) = 


u + (x,y) 


u (x,y) 


(V«» y ■ 0 

J i 

x > 0 
x < 0, 


d 

where by definition f n (u,y) = £ f (u)v.(y), with v(y) normal to S. 

U i=0 1 1 

The variables t and x will play the major role in our expansion 
with y^ as parameters varying in the compact set |y| Rq for some Rq. 
The first term in the expansion will be given by the solution to the Riemann 


problem 
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(1.3) 


u t + ( f 0 (u,y) ) X = 0 

{ u + (0,y) , x > 0 

u ( 0 , y ) , x < 0. 


If the system is strictly hyperbolic and the initial jump is small, the 
solution to the Riemann problem, due to P. D. Lax, is given in [1]. His proof 
involves the construction of the map U(y , , • • • , e^) : Ip if 1 , with y as a 

parameter, U(y,0,***,0) = u_(0,y). U(y , , • • • , e ) represents the state 
obtained by starting from u_(0,y) and travelling z time increments along 
the appropriate shock, rarefaction, or contact curves. Lax obtains the 
solution by showing that U is invertible near e = 0. The solution u can 
be expressed as u(t,x,y) = h(py) with h(X 1 (u_),y) = u_ and 
h ( A n ( u +)>y) = u + * The result immediately extends to the case with multiple 
eigenvalues in linearly degenerate fields if there is a basis of eigenvectors. 

Our result in this paper is 

Theorem 1. Given u_(x,y), there exists e* > 0 small and C > 0 

large, depending only on u_, f ± , such that if u + (0,y) = U(y, , • • • , e n ) , 
U(y,0) = u_(0,y) satisfies 

(a) I e ± (y) | £ e*, 1 = 1,2, “sn 

(b) _if_ p is a genuinely nonlinear field then 
either 

(1) e p (y) * 0, for |y| < R Q 
or 
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d") e p(y) = o, for I y I < R q , and 
Up( u _( 0 *y))*( u _) x <0»y) - ^ p (u + (0,y))*(u + ) x (0,y)| > Ce*, 


then we can construct a convergent power series which is a distribution 
solution to (1*2), 


The solution consists of regions of analyticity separated by shock, 
contact, and rarefaction waves corresponding to the ones in the Riemann 
problem as well as sound waves corresponding to shocks of zero strength in the 
Riemann problem (the case =0). It therefore gives a precise description 
of the singularities propagating from the initial discontinuity (see Figure 
1 . 1 ). 


2nd field 
rarefaction 



1st field 
shock 


2nd field 
rarefaction 


3rd field 
sound wave 



/ 


/ 


/ 


/ 


full problem 


3x3 system 


Figure 1.1 
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Condition (b) prevents shocks or rarefactions in the Riemann problem from 
degenerating to waves of zero strength within the parameter domain |y| < R Q , 
unless they are identically of zero strength. The difficulty with transitions 
to sound waves is caused by the fact that the two flat characteristic surfaces 
joining together in the Riemann problem will not necessarily ensure that the 
two curved characteristic surfaces in the full problem will likewise overlap 
one another. 

One can distinguish between two types of regions, the ones in the 'gaps' 
between waves where the solution is analytic in x and t and the ones in 
the rarefactions where it is analytic in the variable x/t. However, unlike 
the rarefactions in the Riemann problem, this last region is not a simple 
wave, in that characteristics are not flat and the solution is not constant 
along them. 

The regions are separated by unknown surfaces where we impose the 
following boundary conditions: At rarefaction and sound surfaces we impose 

continuity across and given the existence of the coefficients of the expansion 
derive that the surfaces are characteristic as formal power series. Here we 
need condition (b) (1') to be able to determine the sound surface coefficients 
uniquely. At shock surfaces we impose the Rankine-Hugoniot conditions. At 
contacts we impose continuity of Riemann invariants and that the contact 
surface is characteristic. If the contact has a multiple eigenvalue, there 
will be less than n equations imposed. Nevertheless, it can be shown that 
they imply all the Rankine-Hugoniot conditions across the contact. 

The problem (1.1) with initial data restricted to ensure the formation of 
only one shock has been previously studied by A. Majda in [5] where the first 
existence result for such systems with discontinuous initial data is given. 
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Theorem 1 answers a conjecture of R. D. Richtmyer on existence of 
solutions to hyperbolic systems of conservation laws with piecewise initial 
data [6]. 

The proof consists of two parts. First, the coefficients are determined 
and estimated and, last, the expansion is shown to converge. In the first 
part we make appropriate changes of variables (Section 2) which in the end 
only amount to rearrangements of power series. One could, just as easily, 
determine the coefficients of the original variables uniquely, but he would 
face enormous difficulties at the convergence step. To obtain the 
coefficients, we must solve algebraic equations in the gaps, (n-1) linear 
ordinary differential equations coupled with one algebraic equation in 
rarefactions and coupling boundary equations at the shock, contact, 
rarefaction, and sound surfaces. This is accomplished in Section 3. To show 
convergence we use the estimates obtained in Section 3 to carry out the 
majorization argument in Section 4. 


2. Expansions 

Differentiating in (1.2) we obtain 


( 2 . 1 ) 


Uj. + A(y,u)u x + B(u)*u y = 0 


with 


X 


P 


i 





Let A have m distinct eigenvalues < X 0 < ••• < X and let 

12 m 

i = have multiplicity y^ and correspond to the linearly 
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degenerate fields. We choose a basis of eigenvectors so that VA *r = 1 in 

the genuinely nonlinear fields and l r p^l = 1 in the linearly degenerate 

fields. If > 1 then there is a choice to be made in picking a basis for 

that eigenspace. We will adopt the following convention: In a linearly 

degenerate field, r_ will refer to any of the eigenvectors r ,,•••, r 

Pi p^l’ p i ,y 1 

m 

that span the eigenspace. Similarly, in the expansion u = £ a. r., 

j=l J 3 

ar = 'i a , r , , and a will refer to any of the components 
P i P i k=l P i’ P i’ p i 


V 1 ’ 




P i ,P i 

Consider a gap (Figure 2.1) bounded on the left and right by 


<J>(t,y) = x.(y)t + l 4> m (y) t m and \Kt,y) = A^CyH + 1 ^ n (y)t n 

^ m=2 ^ n=2 

respectively. We change variables as follows: 


x = <K5,y) + f(n,y) 


t = 5 + n 


y = y. 


where 5 , a > y are the new gap variables. 

As shown in Figure 2.1 x = <|>, x = are mapped into a = 0, £ = 0 
respectively. 
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The end gaps, the first and the (m+l) st (Figure 2.2) are bounded by only 
one unknown surface. 



Figure 2.2 


Let X A be a fixed number depending on u_ and f^ , i = 0, • • • ,d . We 
will later specify how large X* is. 

In the first gap we let 

x = <Kh,y) + X A £ 

t = n 

y = y. 


Similarly in the (m+l) st gap 
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x = iKS,y) + X*n 

t = 5 

y = y- 

We obtain from (2.1) 

(2.3) u + “r“ (A - <{> - <fi *b)u- + B*u = 0 

n X* v Y n Y y J £ y 

and 

(2.4) u c + (A - i|>_ - t|> *b)u + B*u = 0. 

K X A ^ 5 y 1 n y 

For a rarefaction bounded on the left and right by $ and ^ 
respectively (Figure 2,3), we change variables as follows: 

X - <Kt,y) 

<f> - * 

t = t 

y = y. 

Remarks : The Riemann solution was an analytic function of x/t in 

rarefactions. Expanding the formula for s above we get 
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s 


x/t - (j>/t 
(<MO/t 


l 



X 

t 



0(t), 


so s behaves very much like x/t. 

The transformation above maps x = <K x = to s = 0, s = 1 

respectively. 


t 



Figure 2.3 


In the new variables, 


(2.5) 


Remarks : 


A = A(u,y). 


(t " <(>)u t + {A - ($ t + s(t - <}>) t ) 

- ('f’y + s(i|> - | j ) )y ) *B}u g + (f - <|>)B»Uy = 0. 
As before u new (s,t,y) = u old (x,t,y). Also, in (2.5) 
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The solutions to (2.2), (2.5) are linked by boundary conditions, 
are four types of boundaries: rarefaction, shock, sound, and contact. 

At a rarefaction surface we impose continuity, 

u old^ (T,y),T,y ) = v old(* (T ’ y) ’ T ’ y )‘ In the new variables we 8 et 


( 2 . 6 ) 


u(0,T,y) = v(0, x,y) 

u(l,x,y) = w(r,0,y) (see Figure 2.4). 


<f> (T ,y) 



Figure 2.4 

At a shock surface we impose the Rankine-Hugoniot conditions, 
and u on the left and right of a shock surface <f>, they are 


There 


For v 


(2.7) 


<l» T (u - v) - (f Q (y,u) - f Q (y ,v) ) + <j> y (f(u) - f(v)) = 0 
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with u = u(x,0,y), v = v(0,x,y). 

At a sound surface we impose continuity. For u and v as above and <f> 
a sound surface, we get 

(2.7a) u(x,0,y) = v(0,x,y). 


Consider now a contact surface <J»(t,y) = A p t + ••• in the linearly 

degenerate field p^. To obtain weak solutions we should impose (2.7); 
however, if p^ > 1, we expect that not all of the n equations in (2.7) are 
independent . 

For each x, y we form the normal flux -f^ + $ f and the 

corresponding map U(e 1 , • • • ,6^) , analytic in e and built from rarefaction 

and contact curves only, such that U(0, •••,()) = u(x,0,y). Then v(0,x,y) is 

connected to u(x,0,y) through a p^ contact if and only if v(0,x,y) = 

U(0,***,e ,0, ••*,()) for some e = fe ,,*•*,£ ]• A Riemann invariant 

P ± P ± »l p i’ M i 


for the 


p i 


th 


field is a function R(u) such that V R*r =0 or 

u P-. 


equivalently r(u( 0, • • • ,0 , e p ,0, •••,())) = constant. There are exactly 
n - p^ independent Riemann invariants. To obtain one such set we consider 
Z^(u) , • •• ,Z n (u) the coordinates of the inverse function of U and let Rj = 
Zj , j * p^ Furthermore, we see that VR^ at u = U(0,0, ••*,()) and 

j .k * Pi* 


Lemma 1: The n - p^ + 1 conditions 

(i) ♦_ " A (u,y , 6 ) 

L * 

(ii) Rj(u,y,<J> y ) = Rj(v,y,(j) y ) , j ^ Pi 


(2.7b) 
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imply that 
(2.7) hold. 


X (u,y,<|> ) = X (v,y ,<f> ) 
p i y p i y 


and that the Ranklne-Hugoniot relations 


t* h 

Remarks : X (u,y,io) is the p^ tn eigenvalue of 


Rj(u,y,w) is the j 


th 


Riemann invariant for the flux 
hence analytic in all of its arguments. 


-A(y,u) + o)»B(u) . 
-f Q (y,u) + <o*f(u) 


and 


Proof: Fix x, y and let v = v(0,x,y). Then v = U( £,,•••,£ ) for 

some , • • • ,6^; consequently, (2.7b) (ii) = 0 for j * p^. This 

means v is connected to u through a contact. Therefore 

X (u,y,<|> ) = X (v,y,<p ) and X (u - v) = (-f (y,u) + <|> f(u)) - 
p^ y y p^ ^ y 

(-fg(y,v) + (j)^ f(v)). The result follows from (2.7b) (i). 


The Euler equations have two well known Riemann invariants for the middle 

2 2 2 

field. They are the pressure p = (y - l)(e - (p(u^ + v + w and the 

A 

normal velocity u = E u + 5 v + £ w with (£,£,£) the normal to the 

x y z x y z 

surface. If, as before, the surface is given by x = <J>(t,y,z) we have the 
following three conditions at the middle contact 


CD ft = ~ U Q + ^ V 0 + 4 Z W Q 


(ii) 


" U 0 + ^y V 0 + *z W 0 = U 1 + ^y V 1 + K W 1 


p(p 


0’ V v 0’ "O’ 




e o ) = p(p l’ 


1’ w l’ 


V 
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with (p Q , Uq 9 Vq , Wq, e Q ), (p^ the left and right states. 

One easily verifies that (i) and (ii) above lead to the Rankine-Hugoniot 
conditions for the Euler equations. A tedious computation gives the 
eigenvectors 


— — 


— — 

1 


1 

u - 5 c 


u + 5 c 

X 


X 

v - 5 c 

, r 0 = 

V + £ c 

y 

* 3 

y 

w - 5 C 


w + £ c 

2 


z 

A 


A 

H - uc 


H 4* uc 

— 


— _ 


corresponding to Aj 
and the sound speed 


c 


u 



= u + c, where the total enthalpy 
. It can now be easily verified that 


H 


e Hh p 
P 


Vu»r^ ^ u * r 3 


Vp*ri Vp'r*: 


is invertible. 

We seek power series solutions of the following form: In gaps 


u(£,n,y) = l u (y)c m n n 


ra,n>0 


mn 


( 2 . 8 ) 
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whereas in rarefactions 

(2.9) u(s,t,y) = l u, (s,y)t k . 

k>0 

The first term in both series is obtained from the Riemann solution h(s,y). 
In (2.8), UQg(y) = h(X < j ) ,y) = h(X^,y) since h is constant in its first 
variable in gaps. In (2.9), ug(s,y) is the solution to 

(A(u 0 ,y) - (X + + s(X^ - X + )))(u 0 ) s = 0 

which is the zero order relation obtained from substituting (2.9) into 
(2.5). It follows that 

u 0 = h ( s <\j, ~ V + V y ^ 0 < s i. 1 

and that if we have a p rarefaction 


( 2 . 10 ) 


V U 0’ y) = \f. + s(X t " V 


( Vs = (x * _ V r P ( v y) 


The various unknown surfaces have expansions of the form 


♦(5,y) = x 4 (y)5 + l My)^ 

9 k>2 K 


<Kn,y) = X (y) n + l a (y)p 
9 k>2 K 


( 2 . 11 ) 
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Substituting (2.8) and (2.11) into (2.2) and collecting terms for £ m n n , we 
obtain the following recursive relations 


" A(u 00 ) l (m + 1)u m+l,n + + A(u 00^ (n + 1)u m,n+l = F mn 


with 


mn 


j^n " X ^ 1 " ( A < u > “ A ^ u qo^ ) + < + + *> y B J u f 


( 2 . 12 ) 


+ r X <f, “ + ( A ^ u ) " A (u nn )) “ + < f’)„ B l u 


00 ' 


y n 


+ (4» n - <(> 5 )B*u y 


mn 


Remarks : F mn contains coefficients of u of order < m+n where by 


definition the order of u mn is m+n* 


m 


If we let u mn - <<- 1 ) ml (y)r l ('J 00 ,y), we obtain 


(2 


.13) <«*1>[\ - + = < J l>„ 


uhere <F 1>«. - V (F) mn' 1 ■ 1 - 2 "— •” and V u 00’ y) > r l (u 00- y) > V u 00’ y) 

are the left and right eigenvectors and eigenvalues of A(u 0 Q,y). For the end 
gaps (2.3) and (2.4) we obtain 

(2.14) (n + l)(o ) .. + (m + 1) -J. [X, - X. ](«.)., = (F, ) 

i m,n+l A L i <J) J i'nrt-l,n N i mn 


with 


<f i).„ 


- - A. 


if | <X 4> " V + ^ A(u) " A(u 00^ " *y B J U ? + 


Bu 


mn 
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(2.15) (m +!)(<*.) 


i'm+l ,n + (n + 1} 


X. - X. 

i <fi 


^ a i^m,n+l ^Vmn 


with 


(F.) mn = -l, 
i mn l 


XT ! ( S " V + ( A(u) " A(u 00^ " *y B | U n + 


Bu 


mn 


Remarks : The reader will note the omission, for simplicity, of an index 

on the L's and F's signifying the gap we"re in. 

To obtain the equations in rarefactions we substitute (2.9) in (2.5) and 
collect the terms involving t^ to get for k > 1 


(2.16) 


(x t “ V kU k + A(u 0 ) " ( \ + s(X * ” v) pVs + A ' (u o )u k (u O ) s 


- Oc+1) (4» k+1 + s(^ k+1 - 4> k+1 )) (u 0 ) s = F k 


with 


F k = 


|(* ~ t) - (<}> - X^ t) u fc + (A(u) - A(u q ) )(u - U Q ) £ 


(2.17) 


+ (A(u) - A(u q ) - A'(u q )(u - u Q ) ) ( u 0 ) g 

" (* t ■ % + s ( ( + t " V " (<f> t " v )(u " v* 


+ ('P ~ <l>)Bu y - (<)) y + s(i|) - <f>) y )jBu £ 


Remark : <\>^ and are t ^ ie highest order coefficients occurring 


in F k . 
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m 

If we substitute u k (s,y) = J (a i ) k (s,y)r i (u () ,y) in (2.16) and 
(2.10), we obtain 


(2.18) 


- V k<a A + (VV ' 


m 


+ ( (X 1 ' 


- < k + D (* k+1 + s(t k+1 


s k 

+ J 1 -(A'(u 0 )r j u 0s )J(«j) k 

- w ( v k - 


To simplify (2.18) we note that by differentiating Ar^ = r^ 
respect to u r , multiplying on the left by on the right by r^ 

summing 


*i( A ' (u 0 )r j u o ) = 


■ X A><\, - 

<P P 


V*i 


Jr r. + (X - X )<5 VX *r . 
P J ’I' <t> ip P J 


with Jr 

P 


3r 

, the Jacobian derivative. Since 


£ 


dr . 

1 

i ds 


V Jr j ( Vs 


- W^jS' 


instead of (2.18), we now have 

(x i - y(\) k + «»♦ - v<°a + a - y i 

(2.19) 

= ^ F i^k’ i ^ p. 


use 


with 

and 
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( 2 . 20 ) 


(k + 1)<X + - X^)(a p ) k + (X^ - X + ) (VX p .r j)( a j)k 


- (k + DCX^ - X + )(+ fcfl + s(4» k+1 - <|» k+1 )) = (FJ 


P k 


where 

Bij = (X ± - X p )A i *((Jr j )r p - (Jr^). 

Remarks : Note that = 0 which is why we let j * p in the sum in 

(2.19), Since doesn^t appear in (2.19) we have a partial decoupling 

which will prove to be helpful. 

We now turn to expansions at boundaries and use equations (2.6), (2.7), 
(2.7a), (2.7b). For rarefactions, from (2.6) 


V°’ y) = v 0k (y) 

u k (1 ’ y) = w k0 (y) - 

For shocks we substitute series for u, v, <j> into (2.7) and collect the 
coefficient of x k to obtain for k 1 

(2.21) (k + lH k+1 (u 00 - v Q0 ) + X^(u k0 - v Qk ) - (A(u 00 )u k0 - A(v 00 )v 0k ) = g k 


with 
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(2.22) g k " " U+ x * X 4*> CCu " u 0Q ) " (v " v 00 » 

- [(f Q (u) - f 0 ( U(J0 ) - A(u oq )(u - u Q0 )) 

- (f Q (v) - f 0 (v 0Q ) - A(v qo )(v - v 00 ))] 

+ <J> y (f(u) - f(v))} k . 

Remarks : The zero order coefficient of x is 

X <f>^ u 00 “ v oo^ ” ( f </ u o<P ” f o (v oo ) ) = 0 

which is the Rankine-Hugoniot condition for the zero order Riemann solution. 
To simplify (2.21) we recall that for a p shock 


= U(e 1 , 


,••*,0) and 

11 

o 

o 

> 

^( E 1 ’ £ 2 ’ * * * ’ G p-1 *° ’ * * * ,0 > With 

u( £ 1 

•’ e p-l’° 

,0, •••,()) = 

r p (v 00 ) 

[ 1 ] . As a result we have 

(1) 

*i* v 00* 

•r.(u 00 ) - 0(e p ), 

if i * j; 

(ii) 


,r i (u oo ) = 1 

+ 0(e p ) 


(Hi) 

V v oo* 

* (u oo " v oo ) 

= e + 
P 

°< e p> 

(iv) 

V v 00> 

* (u oo “ v 00 ) 

" °< £ p> 

, if i * p. 


\o ■ £ <VkO r j (u 00 ) 

v 0k " | r J (v 00 ) 


By substituting 
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in (2.21) and using (i), (ii), (iii), and (iv) above, we get 

< k + 1H k + i 0(£ p> + (S - V u oo ) ) ( “iW 1 + 0< V ) 

tot V W’WXVko 

• l x * - YW^Vok ■ <*iV 1 * P 

and 

<k + IH k+i( £ P + o(e p>) + ( x + - x P < u oo ) ) <o p ) ko(‘ + o( y) 

+ 0 ( V Jp (X * ' X j (U 00»<»J>kO 
■ ( x + - y'WHv ■ <* P v 

Since 

x p (u oo > " x p (v oo ) + 0(e p > 

V V“oo> = -^ + °< £ p> 

vvv-V'i^ tn 

we get 

< k * 1)V, 0(Pp) 1 (% - - ( 6l ) 0k ) 

+ %> Y<°>kO - 


i * P 
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(k+ 0 <O) + 


p' p ' p'kO 


+ 0(e P )T * ( Vok " Vk 

where S. *» (S . , ,S.„ , • • • ,S. ), j = T = (T.,‘»»,T ) are vectors bounded 

3 Ji J m i m 

independent of near zero. S and T will change in the next equations, 

but they will remain bounded. Solve for (k + IM^+I in the second equation 
to get 


(2.23) 


(k+ 1H k + l ■ s p- (a \o + T -< 6 >0k + 


1 


2 

e + 0(0 P 

P P 


Substitute in the first equation and divide by A - ^( v qq) t0 get 


(2.24) (« i)k0 - ( 6l ) 0k + 0(0 p )S 1 -(«) k0 + 0(0 )T.(« 0k - P 1 -(g) k , i * p 


where are bounded independent of 

Let us now consider a sound surface 

x * = x p (u oo ) = X p (v oo ) and x and * 
2.5). Expanding (2.7a) we simply get 


as well. 

(j> = X^t + , • • • , where 

are neighboring surfaces (see Fig. 


u k0 (y) = v 0k (y) 


which in coordinates, since uqq = Vqq, gives 


^VkO ( Vok' 


(2.24a) 
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Figure 2.5 


The surface coefficients <j> can be recovered from the gap relations 

(2.12) . Focusing on the equation, if we first let m = n = 0 in 

(2.12) , (2.13) we get 


<x * ' V ( Vi0 ■ ' <x * • VV B(u oo )(u ooV 


Similarly, in the left gap 


(X - X )(3 ) A1 = -( X , - X )l »B(v nr . )(v nn ) . 
P x P 01 4> X P 00 00 y 


Since UQg(y) = v 00^^ ^°^ ows that (a^)^ = ^p^01 satis ^ ie< ^ as a 

result of the gap equations. Next, letting m = k-1 > 0, n = 0, note that 
from (2.12) 

F k-1,0 “ k +k U 01 + k +k B 0 u 00y + F k-1,0 


where F n contains only lower order coefficients of <(>. Hence, from 

kC i j u 


(2.13) 
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k<X * - V'VkO = k *k V <U 01 + B 0 “ 0 0y> + »p ? k-l,0- 


Referring to the original variables, x and t, for a moment 


V u oi = vK (0 * 0) + \ u x (°>°)] 


= £ D*t _A(u 00 )u x (0 ’ 0) + - B Kr>nnJ 


fx' ’ 7 v 00 7 “00y- 


= (A. - X )Z »u - Z *B. u AA 
^ p p x P 0 OOy 


where we denote u = u (0,0) in the gap. Hence, 


(2.24b) 


(X t ' V (k( \\0 - k *k *p-"x) ' ? k-l,0' 


Similarly, in the left gap 


( 2 . 24c) 


<V 


x x )(k(B p ) 


ni - kcf>. Z *v ) « F 
Ok Y k p x J 


o,k-r 


Therefore if Z p (u x - v x ) * 0 the boundary condition (a p ) kQ = (3 p ) 0k will 

determine <j> . To establish Z (u AA )(u - v ) * 0 we let A be a shock 
K p Uu x x T 

surface, for example, and w the function in the gap to the right of u. 
Expanding the shock relations 


<f > t (w(d»,t) - u(<|>,t)) - [f Q ] 


+ d> [g] . 

jump y jump 
at <f> at <|> 


= 0 


and collecting first order terms we easily obtain 
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i (u nr .)u = l (w-.-.Jw + 0(e.). 

p 00 x p 00 x * 

Crossing a rarefaction will yield the same estimate by switching the sides and 
therefore reducing it to the shock case. Crossing a sound surface 
u(<f>,t) = v(<f>,t), again gives the estimate above and, therefore, going through 
all the boundaries 


V U 00 )u x = * p ( u + (0>0))(u + ) x + 0<e*> 

A p (v 00 )v x = £ p ( u - (0 »°>) (u ->x + 

If C is large enough in (b)(l') of Theorem 1, we obtain the desired 
condition. 

It remains to expand at contacts in (2.7b). Equation (2.7b) (ii) yields 


7 u R j (u 00- y ' 0) ' u k0 - V v V v 00- y > 0)v 0k ■ (L j ) k 


where 


( L -j) k = - (Rj(u,y,<t> y ) - R j( u 00 »y,o) 


(2.25a) 


- VR^ (u Q0 ,y»0)(u - u Q0 )) - (Rj(v,y,* y ) 

- R.(v 00 ,y,0) - VR.(v 00 ,y,0)(v - v Q0 ))J fc . 


Here we used the fact that Rj(uQQ,y,0) = Rj(vQQ,y,0), i.e 
are connected through a contact. If we let u^q = £ ^ a i^ k o 


, the initial data 


r i (u 00>» 
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V 0k = ^ ^i^Ok r i^ v OO^ and USe formulas U)* (ii)* (iii)* and (iv) derived 
for the shock expansions we get 


jp. 7R j (u oo > ' t i <u oo ) ( ( “i ) ko • (f, iW + °( e Pl ) s j (w 0k ■ <L j>k 


with S. bounded independent of e ->-0. Since ( VR. (u rtrt ) »r . (u~„ ) ) , 
J P i v j 00 l 00 


i, j * is invertible 


(2.25b) 


( “l>k0 - <e i>0k + °( e pJ S i (W 0k ■ V (L) 1 


with S_p bounded matrices and 


(L) k = ^ (L j } k’ j * p il 

From (2.7b) (i) 


(2.25c) (k + l)* k+1 - ^ (V u » Pi <«00>y' 0) ' r j <u 00 > ) <o j)k 0 * ( L pJ k 

( L Pi)k ' { X p t (tt - y >V - X Pl (u 00' y - 0) “ 7 u X p 1 (u 00- ), ’ 0) ’ (u - u 00 ) } k 


for k > 1. 

At this point we should be able to show that all coefficients can be 
uniquely determined from the formulas established so far. We will do it in 
the next section. To conclude this section, we derive from the conditions 
already imposed that the rarefaction and sound surfaces satisfy a 
characteristic equation. 
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We have 

Lemma 2: Suppose there is a unique formal power series solution. If <(> 

is either a p-rarefaction surface or a p-sound surface, then 

(2.26) (<t> t (t,y)) k = (* p (u,y,$ y )) k = (* p (v,y ,<f> y ) ) k 

with u, v the solutions near <j>. 

Proof : We give the argument for rarefactions, the one for sound surfaces 

following the same lines. 

Suppose u is the function in the rarefaction to the right of <j> and 
v is in the gap on the left. 

Let 

u old <x *t.y> = 

v old (x,t,y) = G(pt,y). 

Then u new (s>t »y ) = H ( s ( * t + and H(o,t,y), G(o,t,y) satisfy 

(2.27a) tH t + (A - aI)H a + tBH y = 0 

(2.27b) tG t + (A - oI)G 0 + tBG y = 0 

and 

H (pt.y) ■ G(pt,y). 
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Differentiating and multiplying by t 


tH t + H„(* t - 1) - tG t - 0 o [* t - i) 


4> 4 

y i y i 

H + H — - = G + G — 
y ± a t y ± at 


Using (2.27a) and (2.27b) with s = the first equation leads to 


H 


,(*t - I-) - Ca - - tBH - G 0 (* t - A) - (A - i)G 0 - tBG 


The second equation, after multiplying by tB^ and adding, yields 


<j> BH + tBH = d> BG + tBG . 
y o y y a y 


Hence we obtain 


(♦ t - A + <|> y B)H a = (4> t - A + <p y B)G a . 

Multiplying on the left by Z p ( H (£,t ,y ) ,y, <f> y ) = ^ p (G(^,t,y),y,<|. y ), we get 

(♦t " V U,y ’V^V (H a " G a ) = °* 

Now (^ p *(H a - G o )) q = ^ p (h(X (J) ,y),y)*h o (X ({) ,y) = 1, where h(a,y) = H^.O.y) 
is the Riemann solution, and since G is in the gap (G o )q = » 0 » y ) = 0* 

Also h a (X^,y) = r p (h(X^,y) ,y) . Therefore, since 

k-1 

(A - X (H,y , 4> )), + T <* — X > (Jt »(H - G )1 = 0 

VT t p v ,7,Y y ' J k vr t p w'- p a (r'k-p 
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and (<{> t - A )q = 0 our claim (2.26) follows by induction. 


Remark : Note that H(y,t,y) = u(0,t,y). 


Expanding (2.26) we get 


+ 2(J>2 t + 3<|>2 t 2 + ••• + (k + 0<|>k + i ^ + *** 


= A (u n ,y,0) + V A • (u - U-. ) + L (u ,y , <}> ) 
p 0 up 0 p ,J ’ T y 


where Lp(uQ,y,0) = 0 and Lp is quadratic in (u - Uq). Hence, for k 1 


we get 


m 


(2.28) (k + lH k+1 (y) = I (V u A (0,y)*r i )(a i ) k (0,y) + (L ) fc . 


i=l 


Remark: (L (u,y,6 ) ), contains u, ,*4), as highest order 

v p y 'k k-1 1 T k 


P y 

coefficients. In fact 


( 2 * 29 ) (L p } k = (y u > y >V - X p (u 00’ y ’ 0) " v u X p (u 00’ y ’ 0) * (u " U 0(A* 


Formulas (2.28), (2.29) are, in fact, expansions valid for rarefactions, 

sound, as well as contact surfaces (see (2.25c)). 


3. Linear Estimates 

In this section we derive a priori estimate for the linear system of 

t* Vi 

equations satisfied by the k u order coefficients with inhomogeneous terms 
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depending on lower order coefficients. These estimates will help in 
determining the coefficients uniquely and subsequently in showing the series 
converges. 

In (2.13), to obtain coefficients of order k, we take m + n + 1 = k 
(k >_ 1). As (2.20), (2.23), (2.25c) suggest, we would expect to determine 
( f > k+1 for shocks, contacts, and rarefactions at the same time we determine 
a^'s. For sound surfaces we can only determine (j> k from the boundary 
conditions, but (2.28) shows that, once determined, the coefficients <f> can 

K, 

be estimated at the previous step. 

Consider the diagram in Figure 3.1 showing the m fields with the gaps 
between them. We let dotted lines signify the various waves. For example, in 
Figure 3.1 we collapse a p-raref action to a dotted line with arrows pointing 
at corresponding faces. 

We now want to consider the coefficients in the gaps at the boundaries of 
the gaps. From (2.8), are the coefficients of the expansion at n = 0, 

the left boundary of the gap. Similarly, a Qk are the coefficients at the 

right boundary. In the first gap we only consider a , at the right boundary 
and in the (m+l) st gap a^, the left boundary. 



Figure 3.1 
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Since there are m boundaries, each with 2n components on both sides, we 

get a total of 2mn unknown boundary components. They satisfy a linear 

system of equations given by the gap equations (2.13), (2.14), (2.15) and by 

boundary equations: rarefactions (2.20), shocks (2.24), sound (2.24a), and 

contacts (2.25b). Note that the p^* 1 equation at each boundary determines 

the surface coefficient. For example, at a sound surface, ^ is determined 
t - v» 

from the p equation at the boundary and, in view of (2.24b), (2.24c), it 
can be solved in terms of lower order terms and hence substituted back into 
the equations (2.12) for the neighboring gaps. (For k = 1 the p t ^ L 

equation is satisfied automatically.) As a second example, the p*^ equation 

at the boundary of a rarefaction region (the continuity condition) determines 
(a p ) k (0,y), (a ) k (l,y). They, in turn, determine $^+1 and \+l 

evaluating (2.20) at s = 0 and s ** 1 . Fortunately, as we mentioned in the 
remark after (2.20), we can solve (2.19) independent of (a ), n . We may 

therefore only consider the n - y^ equations at the p boundary. If p is 
a genuinely nonlinear field y = 1 and we have n - 1 equations. The total 

p 

number of equations for the 2mn unknowns is 


m 

(I n - y ) + (m + l)n = 2mn. 
i=l 1 


The first terra above gives the total number of equations from boundaries, the 
second one gives the total number from gaps. To show the system has a unique 
solution it suffices to prove the linear mapping is one to one. This will 
follow from the estimates ahead. 

We now divide the unknowns into two groups a^ and t^. If we are at 


f* Vi 

the p 1 " boundary (dotted line), we count (a^, •••,<* in the gap on 


the 
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left and (a j , • • • , 01 ^) in the gap on the right as part of a and a 
consists of exactly these components. The rest forms S (see Figure 3.2). 



Figure 3.2 

It follows that a has n*(m - 1) components and & has n*(m + 1) 

components • 

We will be able to estimate a from the boundary equations and £ from 
the gap equations. In rarefactions is the characteristic component 

satisfying the algebraic relation in (2.20). Note that since the boundary 
values of belong to $ and not a, they will be estimated from the gap 

equations and not from (2.20). 

A. Estimates from gaps 

We consider the gap between the and (p+l)^ fields bounded by 

surfaces <|> and tJ ; on the left and right respectively. Lax^s entropy 
conditions give A^^Uqq) 2l ^ — ^p^ U 00^* We £ et ec l ua l ities at the ends 

for sound or rarefaction surfaces. From (2.13) we get 
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X, - X 


( V«H,n - X^Vm.n+l + iTT A - ^i^.n 1 - P 


<a i ) m,n+i - (rfrKr^X'iUj,,, + ftt t-^t - (F i>m,„ 1 > p +1 

l © I©’ 


Let 


X i 

V X i 


p i = 


\. 

1 


- X, 


X. 

l 


- x, 


i < p 


i > p+1 


Then there exits p independent of e A near zero but depending on u_(u), 
such that 0 < p^(u) _< p < 1 for |y^| Rq where u + (y) = 

U( e i » e 2 » * * * » e m ) an ^ | e ^| 5. e *» e * to be chosen. We claim we have 


(a. ) 


i m+1 


m+1 /m + n + 1 \ , % 

,n ' p i ( n ^ a i^ 


Ok 


+ p m . ( n + 0(n + 2) ♦ « *(n + m) 1 . 

P i (m + 1) *m* ••2*1 X - X. ' i ; 0,m+n 

i() i * 


(3 i a ) + m-l .(n + 1) * • • (n + m - 1) 1 , v . 

i (m + l)*m*»»2 X^ - X^ * l,m+n-l 


1 1 


] A (t 1 ) 1 f \ j_ A 1 / V \ 

P i (m + 1 ) *m A^ - K i'm-ljn+l X^ - X 1 m+1 ^i'm.n 


for i < p 
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and 


, N n+ 1 / m + n + 1 \ / N 

’ p i ( m ) ( “i> 


kO 


(3.1b) 


n + 1 ) (m + 2)***(m + n) 1 . . 

P i (n + 1) *••2*1 X. - X. l ‘ i? i'm+n,( 

i <)> 

+ • • • + i i ( p \ 

X. - X (n + 1) i'm,n 
i 4> 


for i > p+1, 


This can be proved inductively on m, say, by substituting formulas (3.1a) for 

(ot.) _i_i into the recursion formula for (cr,) 

l m,n+l ' i m+1 ,n 

For notational convenience we let ( F^) be the sums involving the 
F^'s on the right-hand side of (3.1a) and (3.1b). Hence (3.1a) and (3.1b) 
can be written as 


(3.1c) 


r \ ra+l/m + n + h, . . 

<a l>.*l,n ‘ "l ( „ ) < “i > 0,k + (F i>m, n *<P 


/ \ n+1 /m + n + 1\, N n 

( a _j) „.i ~ Pj [ n + (F.) i > p+1 

l ra,n+l l v m J 1 k,0 1 ra,n — v 


which gives 


(3.2) (“i>k,0 p i (o i>0,k + ^i^k-1,0 


i < p 


(3.3) ( a i> 0 ,k p i^ a i*k,0 + ^i^O.k-l 


i > p+1 . 


In the end gaps we have 
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lSt g a P u new (0,?,y) “ u old ( X * ^ ,0 ,y) = 

(m+l) st gap u new (n,0 ’ y) = u old ( X * T1, ° ,y) = u +^ X * n »° • 

Let ( u ± ( x *?*°»y)) k = I ( a ± ) r i (u ± oo ) with u ±00 = u ± (0,0,y). Then from 

' JL lc 

(2.14) and (2.15) we get the following relations for the end gaps, 1 st gap, 
and (m + l) st gap respectively: 


(3.4) 


, \ n+1 

^ a i^m,n+l P i 


( m+ntl )/« ) +<F,> all 

” \-llk 1 "- n 


(3.5) (« ( ) „ - p" +1 (" + " + 1 ) (a \ + (F. ) all i. 

i m+l,n i v n J y + J. i m,n 

X d» ' X i 

Here y ^ and it follows that if we pick X A large enough, 

depending only on u_(y), f^, we have 0< | | _< P < 1 • Hence it follows 

that 

(3.6) 


(a i*0,k p i(“- 1 ) k + ^i^O.k-l 


1 St n 

1 Gap, 


< 3 ’ 7 > '“l’k.O ■ °i(“+ 1 ) k + <F l > lc-l ,( 


(m + l) St Gap. 


Let denote the vector containing all F^'s from all gaps (3.2), (3.3), 

(3.6), (3.7). We note that the components on the left-hand side of (3.2), 
(3.3), (3.6), (3.7) together form the whole of and the ones on the right- 
hand side of (3.2), (3.3) next to the p^ form the whole of a^. 

Hence (3.2), (3.3), (3.6), (3.7) give us 


(3.8) 


|S k l < p k |i k l * l«v> k l + l<--> k l + Hk-il 
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where |$ | = max{|b , |[ denotes the max norm. 




Remarks : = max|F^'s| and by F^'s we understand (^) k g or 

( ^ 1 ) 0 ^ as t ' ie case may be. 

B. Estimates from rarefactions 

We consider (2.19) with 0 < s < 1. Let 


. , ^ V y> ' V y) 

i ( °' y) ' X„(o, y> - M°,y> 


i < p-1 


A. " ** 

A i ( “ >y) * x - > 

i P 


i > p+1 , 


Then if e = e (y) is such that u.(y) = U(y, £.,••*,£ ,•••,£ ) we have 
p P + J 1* p in 

e p(y) = x ^(y) - ^(y) and hence 


(3.8a) 


0 ^ e o — — C o e *’ 


y ± | < R 


with £q = inf{A i ,|y i | <_ Rq} > 0 by (b) (1) of Theorem 1, and Cq dependent 


k / A j ,(a,y)d0 


on u_, fj only. For i < p-1 we use e 


as an integrating 


k / A. ( a,y)do 
0 A 

factor. For i > p+1 we use e . We obtain: 


d_ 

ds 


k / A t 


( °i ) k e 


k / h ± 


k / A ± 


A i e 


B i (a) k 


1 e s ( Vk, 


A. - A 
i P 


for i < p-1 
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d , . 0 

dS - (a i ) k 6 


k / A k / A 


k / A. 


+ A^ e 


1 0 " , 
17 -T' e <F i>k 

1 P 


for i > p+1, 


Integrating, we get 


(ai) k ( s » y ) = («i) k <l.y)*e 


-k / A ± j -k / A ± 


A i V\ ds ' 


i - k / A i 
♦ / e 8 


T-^X- < F i>k ds ' 
P i 


i < p-1 


C^ys.y) 


• k l A 1 s - k ', A i 


C“4) k (0,y)e 0 + / e s 


A i B i‘\ ds ' 


8 ~ k! A , 

W e s ‘ 


X-^TT <F i>k ds ' 
i P 


i > p+1. 


It follows from (3.9) and (3.8a) that 
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l (a iV s,y) l - l (a i>k (1,y) l + c o e * sup 

0<S<1 


+ sup I (F ) I i < p-1 

0<s< i 

(3.10) 

l (a iV s,y) l - l (a iV 0,y) l + C 0 G * sup l B i ,0t kl 

0<s<l 


+ sup | (F . ). | i P+1 

0<s<l 


where 


i - k/ a i , 

(^i>k = ^ 6 X r ~TT * F i*k ds 1 - P_1 

s pi 


s _k / A i 

<F A ■ I c s ' i-hr <F A ds ' 1 > » +l 

0 i p 


and since 8 ip = 0 by (2.20), = (( , • • • .(a^ ) fc • * * ,( 

Now the boundary condition (2.6) gives us (see Figure 3.3) 


P 



Figure 3.3 
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(« i ) k ( 1 , y ) = ( y i) k o (y) 1 - P_1 

(a iV 0,y) = ^i^ok (y) 1 - p+1 

with 

W mn * l <Y i > mn r i < »00 ) 

V ■ l r i (v 00 ) 

where we have 

w oo = u o (1 ’ y) 
v oo = u 0 (0 ’ y) 

by the continuity of the Riemann solution in rarefactions. Now 

\ 0 ( y ) i < p-l, C Qk (y) 1 - P+1 

belong to our vector 

We will now adopt the convention that Cq will denote a constant 
depending only on u_, f^, but it will get larger from equation to equation. 
With this, from (3.10) we obtain 

|(°i) k ( s ’ y >l < |S k l + c 0 e * SU P l\l + SU P l F ikl 1 * P* 

s s 

For £* < we get by absorption 
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sup loj < 21^1 + 2 sup 1(^)1 
s s 

where as before = ^ik’ F 2k’ * * *^ * Usin 6 this in (3*10) we obtain 

|(a i )k ( 0 >y ) i < |* k | + C Q e*(2|^| + 2 sup Kf^l) + sup ) \ i < p-1 , 

s s 

|(a i ) k (l,y)| < |S k | + C Q e*(2|S k | + 2 sup 1(^)1 ) + sup |(F )J i > p+1. 

s s 

Now, by the boundary conditions 

(a i) k ( 0 >y) = ( ^ i ) ok 1 < p - 1 

(a i ) k (1,y) = ( VkO 1 > P +1 » 

with ?o,k’ * — P_1 and \ O’ 1 — P+1 belonging to a k * In fact, counting 
all rarefactions, they are the part of ^ on the faces of rarefaction 
boundaries* We call them a • Hence we get 

O.n) |(a R ) k | < (1 + 2C 0 e A ) |? k | + C 0 sup IfJ. 

s 


C. Estimates from shock, sound, and contact boundaries 

The a k 's occurring on the faces of shocks a^, sound a<, D , and contact 
surfaces a c are handled by (2.24), (2.24a), (2.25b). From (2.24) 
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which, for small depending on u_, implies that 

l^ a SK^kJ - ^ + l^ b SK^kl + C 0^ g ^J* 

Similarly from (2.25b) 

l<a c ) k | < [1 + 0(e*)) |(ff c ) k | + C 0 |(L) k | 

and from (2.24a) 

I ^ a SD^J = I ^ b SD^kl * 

This together with (3.11) gives us 

(3.12) |a k | < (1 + C Q e*) |b k | + C 0 (|(g) k | + | ( L) k | + sup IfJ). 

s 

Combine this with (3.8) to get 

|a k | < (1 + C Q e*)p k |l k | + C 0 (|(a + ) k | + |(a_) k | 

+ Ifk-il + Ks) k l + U L > k l + su p Iffcl)* 

s 

By choosing e* smaller, but depending only on u_, f^ we can make 

(1 + C Q e*)p < V 2 • Hence the above together with (3.8) yield our main linear 

estimate 



(3.13) 
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which holds for our choice of and for Cq depending on u_, f^ only. 

Now a^, satisfy a linear system 



where H k comprises of all the inhomogeneous terms (a + ) k , (cO^, , (L)^, 

g^> The estimate (3.13) shows that the 2mn x 2ran matrix A is 

invertible. Formulas (3.1), (3.4), and (3.5) will give directly the rest of 

the coefficients in the gaps. Given the initial values (a^)^(0,y) we can 

solve the O.D.E. (2.19) for 0^s^l,i*p. We can finally recover the 

rest of the unknowns, a in rarefactions and the surface coefficients, from 

P 

the p^* 1 equation at each boundary. 
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4. Convergence 

In this section we prove the convergence of the power series constructed 
in the previous sections by employing a variant of the technique of 
majorization. To carry out this process we must consider our variables s, 
y complex with 

yen y = {yi€$. d(y i ,[-R 0 ,R 0 ] ) < «, i=i , — ,d} 
s en s = {se$, d(s,[o,i]) < 6}. 

Remarks : 6 is a small number less than 1 to be chosen later and 

d(s,[0,l]) represents the distance from s to [0,1]* 

If we begin with complex analytic initial data u ± and complex analytic 
coefficients in our equation (1.2), it is clear that all our equations will 
hold for y and x complex. 

A, Auxiliary Lemmas 

We define 

\ = {u(s ,y) analytic in n x o sup (d(s,n C )»d(y,fi C )) k |u(s,y)| < “>} . 

y sen s y 

s 

yen 

y 

It follows that are Banach spaces with norm 

|u| = sup (d(s,n^)*d(s,^)) k |u(s,y)|. 

\ sen s y 

s 

yen 

y 
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We will use the notation d g = d(s,fi g ), d y = d(y,Qp and note that d g , 

d < 1 if 6 1 . Hence | u | „ < |u| . 

y k+1 \ 

Lemma 1: (Hormander [3] , p. 117) 


(4.1) 


for u H^. 


e(k + 1) 

\+l 





H. 


k+1 


< e(k + 1) 


V 


Proof: It suffices to consider u(s), sefi and show the first 

s 

inequality. Fix s€ft and let e < d . Then Cauchy^s inequality gives 

s s 

|u'(s) I < e sup |u(?)| < e ^(d - e) k |u| u Choosing e = d / (k + 1) 

| ? -s|<e S V s 

we obtain 


u'(s)| < (k + 1)(1 + k _1 ) k d s k 1 |u|^ < (k + l)ed s k 1 |u| . 

k+1 

The lemma results by multiplying through by d g and taking sup over 

sen. 

s 

Lemma 2 Let C > 0. Then there exists 6^ = <S A (C) such that 
e rC < ^ f for 0 < r < e, 6 < 6*. 
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Proof : Let f(r) = e rC - f . Then f(0) = 0 and 

e , f N „ rC 6 

f (r) = Ce T . 

(6 - r) Z 

We have 

Ce rC < Ce 6C < ^ if 6 < = 6*(C), 

1 6 

and since j < j , f'(r) <0 for 0 < r < 6, 6 < 6*(C). Hence 

(6 - r) 

f (r) <0 for 0 < r < 6. 


Lemma 3: Given N > 0, there exists 6^(N) such that 

<$* - r 

*n C §-- - p ) < -N(r - p), 0 < p < r < 6*. 

Proof : Let x = r - p and e = 6* - r. It suffices to show 

^ lVx/J < _Nx 0 < x < 6*. 0 < e < 6*, 

or x -1 £n( 1 + x/e) > N. 

There exists Sq(N) such that 

x" 1 An(l + (N + l)x) > N, 0 < x < S Q . 

Take 6* = min(6Q,l/(N + 1)). Then x 1 in(l + x/e) > x -1 £n(l + (N + l)x) _> N 
1/e > 1/ 6* >. (N + 1) and x < 6* _< 6^. 


since 
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Lemma 4 : Def ine 


j -k / A i (a,y)da 

/ e S u(s',y)ds'', i < p-1 


u = < 


s -k / A jL (a,y)da 

II e S ' 

0 


u(s" .yjds' , i > p+1 


for u € H. . , k > 1, 
k-1 — 


Then there exists 6* , depending only on u_ , f , such that 


(4.2) 

(4.3) 

for 6 <_ 5^, with Cq, e 



as in (3.8a). 


Proof: (3.8a) will 

hold 

for s ,y € J2 x Q if 

J s y 

6 is small 

depending 

only 

u_, f ^ . It suffices 

to 

consider the i > p+1 

case only. 

Fix sen 

s 

and 

s be the point 

on 

[0,1] closest to s. 

Let r = | 

* 1 

s - s . 

(See 


Figure 4.1.) 



Figure 4.1 
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Then if s' € [0,s ] 


(4.3a) 


Hence 


s 


-k / A 


s 

e 

= 


-k I A t 


-k I A 

JL J- 


-k(s -s')e Q krC 0 e* 
< e e 


by (3.8a), 




s* ~ k i A l 

/ e S 
0 


u(s',y)ds' 


< e 


krC o e * ' Vl 


* > * 
s -k(s -s')e 

k rr • J « 0 d.' 

(fid r 1 0 


rC o e * 

Tk^I 


rC_ e \ k-1 . . 

e 0 * 1 / , -ke 0 s 

1 - e 


ke 


0 


H 


k-r 


Applying Lemma 2, for 6 < 6 ^(Cq we get 


e ^ C ° £ * ill _ke 0 

1 (S - r )*- 1 I I'V " 1 


since s < 1 , 


* 10 * i0 

Next , let s « s + re and s" = s + pe • Then 


X 2 = 


s 

/ e 
* 
s 


-k 


u(s',y)ds' 


r k(r-p)C 0 £* 

< J e 


u 


1 H \-l 

dp* 


fK .k-1 , k-1 

(5 - P) d y 


rC. 


E * |u| \-i , r <*- l > 


f e 

(5 - r) k_1 d k-1 0 

y 


[( r -p) c 0 e* + 


dp. 
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If we let N = 1 + Cq e A , the integral above equals 

, r -(k-l)Cr-p) (k-1) [(r-p)N + ta(-f=-|)] 

J e e 

0 


Using Lemma 3, for 6 < $*(N) we get 


< /V (k - 1)(r - p) dp 

0 

<f (1 - e‘ kr > 



e — ( k— 1 ) 1 3 > k > 1 
, k - 1 

for k > 1 


since r 1 . Hence, for 6 < 6^ with 6^ depending on u_. 

Mil 

|t u(s)| < I. + I_ < 

” z “ (d d 

s y 

The inequality (4.3) follows immediately. The inequality (4. 

mm 1 _y P 

observing that x (1 - e ) < C for x, C > 0. 


-1 


vk-1 


oC. 


ke. 


(1 - 


-ke„ 


Vi 


B. The majorant 

Ultimately, we want to show that 


|<*i>k 


i 


H k _ x (k + 1) 


_ a 

2 k 


(4.4) 


(o ) < - l m + 

v Vmn L - f . 1N 2 , x n 2 ^ n J I 

H k— 1 + ^ ( n + n 

l (k + 1) * k+i lH k . 1 1 ^T7V‘ k 


for i = m, m + n = k, k > 1. 


f^ only 

(1 -e- kS )l. 


2) follows by 
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Here <J> denotes any boundary surface and 

a(z) = Oj z + d- 2 z" + •••, a > 0, 

will be a convergent power series. We are not yet ready to say 
what &(z) is. The & k 's will satisfy a recursive relation which will be 
determined during the course of majorization. 

Remark: We define |(a_,) | = sup (fid ) R ^ |(a,) (y)|, where 

I 1 D * n *H. , yen y I 1 ” I 

k-1 J y 

(a.) (y) is regarded is a function of s and y with s = 0. 

I mn 

Let a Q > 0 be an upper bound for all zero order coefficients 
|< a i>o (s ’ y) |’ | (a i ) 00 (y) |’ I ^I (y) I = l x + < y )|» as well as their derivatives 

j WoMVVyH Wool’ll +l| with s ’ yen s x V 

We now state the Main Majorization Lerama, which will be proved in C, and 
use it for the rest of B. 


Main Majorization Lemma: Suppose we have a o » ^ > # * * ,a k-l ’ ^ 2. ^ so that 

(4,4) is satisfied up to and including the index k-1. Then 


< F i>k 


< — (Q x ( 0 .(z) ) + zR^aCz) ,z) ) k 


(M 


H k _ 1 (k + 1) 

I . (m + n + 1) ! 1 


1 


(4.5) 


i'mn - m! n! , , ..2 , , ,.2 "*2 ' ““2'k 

H k-1 (m +1) (n + 1) 


■ lk h a j. ir J JK 


(Q. + zR ) , m+n = k-1 


H k _ L (k + 1) 

(L ) I < l 

1 K 'H k _ 1 (k + 1) 


2 (( ^4 + zR 4^k 
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with (F i ) k , (F i ) nn , (g i ) k , (L ± ) k from (2.17), (2.12), (2.22), (2.25a) and 
(2.29) respectively . 

In (4.5) Qj(<i) = ^j2 + ^j3 ^ + ***> 3 = 1>2>3,4 is a convergent 
power series beginning with quadratic terms and Rj(a,z) is analytic at 
(0,0) bearing no relation to the Riemann invariants introduced before. 


Remark : Qj , Rj will, as expected, involve majorants of the coefficients 

A,B,L,*** of the original equation. We note that (Q^ + zR^) + (Q^ + zR^) = 
(Qj + Q 2 ) + z(R^ + R 2 ). Hence sums of functions of this form have the same 
form, and we will simply denote them all by Q + zR in spite of the fact that 
they may differ from equation to equation. 


Assuming the lemma we now prove (4.4) for the index k. 

We consider ( ^i) mn from (3.1a), (3.1b). For i p, by (4.5) 


0 j (n + l)(n + 2) » * * (n + j) 1 .. . . 

i (m + l)»m»**(m - j + 1) | | ' i'm-j ,n+j | 

/ J c (n + 1) * » « (n + j) # (m + n + 1)! 

— x 0 (m + l)*»*(m - j + 1) (m - j)!(n + j)! 


1 1 

(m - j + l) 2 (n + j + l) 2 


(Q + zR) k 


, /-m + n + 1 
'0 1 - n 


1 1 

(m + l) 2 (n + l) 2 


pj (Q + zR ) k 


. p j (m + l) 2 (n + l) 2 

i 2 (m - j + l) 2 (n + j + l) 2 


where 


p *p and p , p < 1. 
1 2 X 1 1 2 
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The expression 


p j (m + l) 2 (n + l) 2 

*2 (m - j + l) 2 (n + j + l) 2 


is bounded independent of 


j since we write it as 


p{ a + D 2 
2 


(m + l) 2 (n + 1 ) 2 

(j + l) 2 (m - j + l) 2 (n + j + l) 2 


and 2 is bounded by 4, by considering j < y and 

(j + 1) (m - j + 1) 

j > y • This S ives 


| (F . ) I < £ T ( m + n + 1 )(Q + Z R) 

1 mn H k-1 ” (m + 1) (“ + 1) n 1 


'(l + P- + Pj + 
1 1 


i n ^ 

... + p J 

1 


(4.6) 


r + : + ‘kq, * 


(m + l) 2 (n + l) 2 


2 k 


for i < p. 


Similarly 


(F. ) 

i mn 


H, 


(4.7) 


k-1 


(m +1) (n + 1) 


/m + n + 1 
^ m 


)(Q + zR) 


for i > p+1. 


m, n , 


Remark: In all of the above m + n + 1 = k 
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We now estimate (F^) k 's * n rarefactions from their formulas given after 
(3.10). Using the first estimate in (4.5) and (4.3) of Lemma 4 


(4.8) 


(F i>k 




H k _ 1 0 (k + 1)‘ 


We are now ready to get an estimate on a^ , from (3.13), 
For k _> 1, since a ± is analytic 


Jc-1 


|(<* ) (y)| < A C k 1 < A l - , L— 

1 * k 1 ° ° (k + l) 2 (6d y ) k_1 

where Aq, C, are appropriately chosen and C < Cj. Hence 

|(a ± ) k | < 2 ^ zR(z) )k 

1 (k + 1 ) l k 

for R(z) = A l z 1 analytic at z = 0. Also (4.8) implies 
i>0 

x k-l ,k-l i c | 1 , . 

5 d sup If, | <7 (Q + z R ), . 

7 °< s < 1 0 (k + ir 

This, together with (4.5), (4.6), (4.7), applied to the right-hand side of 
(3.13) implies 


(4.9) 

l (b Vk| H . 

k-1 

with C 

depending on u_, f. 


Vk 


H k _! (k + 1) 


j (Q + zR ) k 



55 


To estimate the remaining coefficients in the gaps we consider (3.1c). We 
get, by using (4.9) 


i ro+1, x 
p i ^ a i Ok 


Vl 


. m+1 (m+ l)^(n+ 1)2 

< P. o — 

(m + n + 1) 


(m +1) (n + 1) 


~2 (Q + zR ) k ‘ 


Letting C, which depends on e^, u_, f^, get larger from equation to 
equation, as we did with Cq, we get 


< 


C 

(m + l)^(n + l)^ 


(Q + zR) k « 


Using (4.6) as well, we obtain from (3.1c) 


^°i*nrt-l,n 


H, 


(4.10) 


k-1 


(m + l)^(n + l)^ 


r + ; + 1 )(q + **>, 


for i < p. 


Similarly 


(a.) 


i m,n+l| — 


(4.11) 


(m + 


l) 2 (n + 


1)‘ 


r + ; + Mcq + z r), 


for i > p+1. 


Remarks : Formula (4.10) holds for the (m+l) st gap too, and (4.11) holds 

for the 1st gap (see (3.4) and (3.5)). 


To get estimates for the rarefaction coefficients we consider (3.9) with 
(s ,y) € fl g x Letting H k _j = x ••• x H^_ 1> (n - 1) times, we 

introduce the map T defined by 
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-k / A ± 

(a i ) k ( 1 ,y)e S + T 1 (A i B^v) + Tj- ~ — (F^), i £ p-1 


(Tv) i = 


P i 


-k J S A. 

(« t ) k CO,y)e 0 + t 1 (a 1 B -v) 4- T ( 1 (F ) ), i>p+l 

P 1 


for v €H n ) . 

k-1 


Remarks: Since in (2.20) 8. = 0, in the formulas above 

i,P 


B. = (8. , • • • , B. , 8. , , , •••,B. )eif 1-1 . 

i 1,1 i,p-l’ i,p+l i ,m 


T .'s are the maps defined in Lemma 4. Note that (3.9) means T (a), = (a) . 

± k k 

We want to show that T is a contraction mapping some ball in to 

rC X 

itself. This will give us a fixed point in the ball. 

Using (4,3a) with s' = 0, and (4.9) we obtain 


(«i)k(°. y )e 


s 

-k / A. 

0 1 


r , rC 0 e * . , rC n e. 

< 2 (Q zR) k* TT — 6 — j e 

(k + ir d k 1 6 

y 


Applying Lemma 2 we obtain 


(4.12) 


-k / A. 


a i (0,y)e 


H 1 

k-1 (k + 1) 


2 (Q + zR > U - 


Similarly 
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(4.13) 


-k / A ± 


“jU.y)© 


H, — 2 ^ + zR ) k * 

Vi (k + i r k 


Since 


| A iV\| H - c o e * Kl 

k - 1 


by (4.2) 


(4.13a) 


| Tj. ( A i B ‘“k^L - C 0 e * 

\-l 


w • 


As always Cq depends only on u_, f^. We choose small enough so that 

C 0 e * < ** Using the first inequality in (4.5) and (4.3) we get 


Vr-hr < F i> k ) 


\-i V k + » 


2 (Q + zR ) k - 


In conclusion, adding all the estimates above 



max 

i 


i <T “k> : 


H 


k-1 


(k + 1) 


2 (Q + zR ^k + C 0 G * 



Choose D so that C + Cq D D, which is possible since C^ < 1, Thus 



< 


D 

(k + l) 2 


(Q + zR) k 
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if 



< 


D 

(k + l ) 2 


(Q + zR) fc . 


Now 


I <TV) 1 ' (Tu, i|„ = |V A i V (v - u) 5|h , < C 0 £ * | v - u | n-l 

Vi 1 Vi 


with C Q e A < 1 as in (4.13a). Hence T is a contraction, which has a fixed 
point ( a ) k satisfying 


(4.14) 


N 

k-l 


< 


D 

(k + l) 2 


(Q + zR) k . 


Remarks : D tends to infinity as e Q tends to zero since C does, which 

means that the radius of convergence of our series approaches zero as 
rarefactions degenerate. 

In (4.14) (a) fc = ((aj)^ * * * > ^ a p_p k > ( “p+pR’ * * * » ( a m ) k ) * As 3 flxed P oint 

of T , Or is the solution to the rarefaction O.D.E.'s (2.19). In Section 3 

we could have obtained the existence of (a) k by solving the initial value 

O.D.E. in the complex domain. However, in this chapter we were able to obtain 

the estimate (4.14) in addition to the existence. 

Although the rarefaction surface coefficients could be obtained from 

(2.20), we cannot prove the desired estimate on them from the equation because 

of the (k + 1) factor in front of the (o ), term. Formula (2.20) is not 

P K 

adequate for bounding derivatives of <f> or ty. Fortunately, we have (2.25), 
(2.28) which were derived as a consequence of <J> being characteristic. 
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From (2.28), using (4.9) and (4.5) we obtain 

(4.15) |(k + 1M, I < j (Q + zR). 

' lH k -l < k+1 > 

which holds for rarefaction and sound surfaces. By (2.25c), using (4.9) and 
(4.5) again, it clearly holds for contact surfaces as well. 

We now go back to (2.20) to obtain the estimate on a p (s,y). We use (4.5) 
to bound (F ) k and (4.14), (4.15) to obtain 


(4.16) 



< 


» ! i> 2 


(Q + zR) fc . 


Remarks : To get (4.16) we needed to estimate 


sup (d d ) 

sen 1 y 

S 


k-l 


+fcfl + s( \+l 


W 


< I Vl + S <Vl - W| < 

sen J 1 


(k + 1) 


yen 


2 (Q + zR) k 


by (4.15) . 


Finally, the shock surfaces coefficients from (2.23) satisfy 


(4.17) 


(k + l)<j> 


k+1 


< o ( Q + zR V 


H k _i (k+D 


Remarks : To get (4.17) we used (4.5) to bound (g ) k and (4.9) to bound 


(a) ko> (e) 0k from ( 2 - 23 > 
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Also, <_ — is incorporated into the constant C. 

e + 0(e) e 0 
P P 

Consider the sum of all C(Q + zR) k 's from (4.10), (4.11), (4.14), 

(4.15), (4.16), and (4.17) and call it Q + zR as agreed. Now set 

(4-18) a k = (Q(a(z)) + zR(<t(z),z)) k , k > 1. 

Remarks : (Q + zR) k in (4.18) contains coefficients of <£ of order less 

than k. By the implicit function theorem the equations a(0) = 0, a = Q(a) + 
zR(a,z) have a unique analytic solution d(z) whose coefficients satisfy 
(4.18). 

With this definition of a k , (4.10), (4.11), (4.14), (4.15), (4.16) imply 
(4.4) for the index k. 

Assuming the Main Majorization Lemma we have thus proved (4.4) for all 
k > 1. 

Remarks ; When k = 1 the Main Majorization Lemma implies (4.5) with 
= 0 and R^^ depending on aQ only. 

C. Proof of the Main Majorization Lemma 

Lemma 5: 


l l - 1 < K ° 

y=0 (m - y + 1)^ (p + l)^ ( m + i) 


with Kq a fixed numerical constant. 
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Proof : We let [x] be the integer part of x. Then 


[?] 

T< 2 l 


m , 

l 1 n _ l 

p=0 (m - p + l) 2 (p + l) 2 — p=0 (m - p + l) 2 (p + l) 2 


(by symmetry) 


00 

1 — 


(m - [j] + l) 2 p=0 (p + l) 2 


°° 1 
l 7 * 


(m + l) 2 p=0 (p + l) 2 


Lemma 6: 




< 



+ m 0 


+ n. 


Proof : Consider + m 2 objects. Then the left-hand side represents 

the number of ways we can choose n^ objects out of the first m^ and n£ 
out of the remaining n^. The right-hand side represents the number of ways 
we can choose n^ + ^ out of m^ + m 2 with no restrictions. Hence, the 
inequality in Lemma 6 becomes evident. 

Lemmas 5 and 6 are among the tools used for a proof of the Cauchy- 


Kovalevsky Theorem in [4] 
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Lemma 7 : Let 


u(£,n,y) = l u (y)S y n v 
y>0 w> 

y>o 

v(5,n,y) = l v (y)5 y n v 
y>0 y> 

v>0 

w(5,n,y) = l w (y)5 y n v 
y>0 y ’ 

y>0 

and suppose that for positive constants , i = 0,»««,2, we have 


0 


1 /y + V\ , . 

T ( v 1 < y+v < m+n 


! UU ’ V IH M+V _ 1 1 '(y + l ) 2 (v+ l ) 2 ' V '~ y+V 


1 


Wy ’ v| Vv-i < + !> 2 < v + 


1 (\1 + V \ / 

7772 ( v ) 6 y+v’ l<M +v <® + n 




with 

*(z) « I a Z 1 a > o 

i>l 1 1 “ 


b(z) = l b z 1 fa, > 0 

i>l 1 “ 


c(z) = l c z 1 c > 0. 

i>0 1 1 
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U 00 * | w oo are also bounded by Tg ag, bp respectively , ag, 

H» H. 


^ | u ool ’ | w oo 
0 

b 0 > 0, then 


(4.19) l(u»w) 


mn 


< — ~ ( m * n )((a(z) + a n )(b(z) + b n ) ) 


Wi (m + 1)2<n + 1)2 “ 


ra+n, 


with m+n > 1. 


(4.19a) In case Uqq = 0 = aQ or Wqq = 0 = bp then we get (4.19) under the 
weaker hypothesis that w , u y V > respectively, satisfy their estimates for 
1 < y+v < m+n-1 only. 


(4.20) 


(u*v) 


mn 


K 0 T 0 T 1 


H 


m+n 


(m + l) 2 (n + l) 2 


[(.(.,♦ Vcw] 


ra+n, 


ra + n > 0. 


(4.20a) In case uqq = 0 = aQ we get (4.20) with ^ u " v ^ mn 


as left- 


H 


m+n-1 


hand side , m + n _> 1, under the weaker hypothesis that ^ satisfies the 


estimate for 0 _< y+v _< rafn-1 only. If, in addition , Vqq = 0 = C ( 


“0 


the 


hypothesis on v could be weakened to 1 _< y+v <C m+n-1. 
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Proof : We have 


( 4 . 21 ) 


|(d y 6) m+n 1 (uw) ]nn | < (d y 6) m+n 1 |u mn | |w 0() | + (d y 6) m+n 1 |u_| |w 


00 mn 


m n 


+ l l KJ (d„ «) 

y=0 v=0 
p+v*0 ,m+n 


. y+v-1 


yv| y 


w I (d 6) n+n-(M+v)-l ( 6) 

m-y,n-v| y ' v y J 


< T a w + T, b. u 

— 0 0 m,n 10 111,11 


ra+n-1 


H 


m+n-1 


- I 


u w 

u v I p,V Ih I m-y ,n-v| 

y+v*0,m+n P+V_1 m+n-(p-v)-l 


T T 

01 /m + n 


(m + ir(n + 1) 


i r ; ")[»„ <> m+n - b 0 «_j 


+ T 


y + Vwm + n - (p + v) \ 1 1 


' T V (P + Vwm + n - ^ii + vh 

01 1 K v ){ n - v J 


y,v 

y+v*0 ,m+n 


(p + l) 2 (m - y + l) 2 


a . b 


(v + l) 2 (n - v + l) 2 p+v ®+n-(y+v)‘ 


The second term above is 
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m+n- 1 ra , - 

<XT V n ]afa f) - - 1 

- 0 A 1 1 n J Z m+n-£ l L n , ^ ..2 , lN 2 J 

1=1 y=0 (y + 1) (m - y + 1) 


( l — — 2 1 2 ) 

v=0 (v + 1) (n - v + 1)^ 


, T 0 T 1 K 0 r m + m+ y 1 , 

1 (m + l) 2 (n + V ) 2 1 n J il 


The estimate (4.19) follows. The result in (4.19a) follows from the fact that 
the right-hand side of (4.21) will not have a w mn or u mn term in case 
uqq, respectively Wqq, is zero. 

In case m = n = 0 (4.20) follows immediately. Hence, we assume 
m + n 1 . Since dy < 1 


(4.22) 


(d 6) m+n (u»v) < (d 6) m+n * u v n _ + (d 6) m+n v u_ n 

y m,n — y m,n 0,0 y m,n 00 


+ M (d ? 5) 

y+v^0 ,m+n 


y+v-i I I (d 6) m+n-(p+v) 

m-y,n-v y 


< u T- C ft + v T. a_ 

| m,n| H 2 0 | m,n| H 0 0 

m+n-1 m+n 


+ 1 1% v| 

y,v 1 V ’ 'H 

y+v^O ,m+n 


m-y ,n-v 


y+v-1 


m+n-( y+v) 


< T 0 T 2 1 (m + n + 1)! r 

- /_ . . \ 2 . , .2 ... .. L 0 i 


(m + 1) (n + 1)' 


- . + c a, 

m+n m+n C 
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+ T„ T, I (P + '•](- + n -_CP + v) + i).^ 


y,v 

y+v*0 ,m+n 


n - v 


- y + 1) 


(p + l) 2 (m - y + l) 2 (v + l) 2 (n - v + l) 2 U+V m+ii-C P— v) 


Since 


r * v )r + - V- v u) + *)<- - - + 1) < r + r l x« ♦ i) - 


(m + n + 1) ! 
m! n! 


(4.20) follows. The result in (4.20a) follows since the right-hand side of 

(4.22) will not have the (d 6) m+n v terra in case a n = 0. 

y ra,n 00 U 

Therefore, (4.22) will be valid with I (d S) ra+n ^(u # v) I as left-hand 
* | y mn| 

side. The rest of (4.20a) is immediate. 


r» y V 

Corollary 1: Let u i = A 11 ’ * = l»***» n satisfy 

hypothesis of u in Lemma 7. Then, for a = ( , • • • , a^) , a multi-index 


(K 0 V 


<«“>. n | < ^ 

’ Vi c " + » <” + » 




Proof : Follows from Lemma 7 (4.19) by induction on |a|. 


Corollary 2: Let u = ( u i>***» u n ) as in Corollary 1. Let 

a(u,y) = £ a a (uQQ,y)(u - Uqq)” be analytic in the variables 
u,y € { |u - u nn | < e} x ft . Suppose, by letting u nn (y) = u(0,0,y) , we have 


l l a J u oo (y) ’ y )l < a «> y€ ^> 


- V 


i > 0 
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(in this case we say a majorizes a). Then 


i(u(5,n,y),y)) mn | H 


c; i 'h* (i o , o iw, L' 


m+n-1 


(m + l) 2 (n + l) 2 


m + n > 1 . 


Proof: By Corollary 1 


l (u - u 00>“U| H 


< T 0 V 


‘m+n-1 


since (u - Uqq)qq = 0. Hence 


(m + l) 2 (n + l) 2 


(" * n H(a(z))l“l) 


m+n 


(d 6 ) 

y 


m+n-l 


< ,i„ . f , | a « Cu 00- 5,) l l (u ‘ "ooC.nL 

’ Z=0 a =£ 1 ’ 'H 

1 1 m+n- 1 


(K 0 V 


s \ uu /m + n\, \- | _ / 

- ^ ~ .2, ~ 772 ^ n ) (a ^m+n ■ l I a (u 00’ 

£=0 (m +1) (n + 1) |a|=£ 


y) 


r; n Xa(K 0 T 0 a)) 


(m + l) 2 (n + l) 2 


m+n 


We are now ready to prove the second inequality of the Main Majorization Lemma 
(see 4.5). 


We consider the terms that enter in from (2.12). By the remark 

following (2.12) and by the hypothesis of the Main Majorization Lemma u, <j> 
and in these terms satisfy (4.4). In what follows ra + n = k- 1^0 

and p < m, V < n. 
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First, we remark that 


yv 


m 

- I 

i=l 


m 


(a i>UV 


r i^ u 00^’ 


V s>y) ' r i <u 0 > 


so that 


yv 1 


< T o 




with T q = T Q (u_,f i ). By (4.4) 



X.) 
v y 


H 


y-l 


(y + l)*. 


l 


y+i 


H 


y-l 


(y + IV 


a 


for 1 < y < m. 


Since i|> doesn't depend on 5, (tf>^ - X^)^ = 0 if v > 0. Therefore 


(4.23a) 


l<* - *,) I < 1 =• ' T C + V )“ 

1 ” V, Vv-l “ <" +1 > (» + 1) V IHV> 


for l<y+v<m+n. 


Similarly 


(4.23b) 


|< X * - Vwv| 


< 


H 


y+v 


(y+1) (v+l)‘ 


ry + vn 

V. V ' 


y+v, 


for l<y+v<m+n 
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Let be the unit vector in R n . Let a be a majorant to 

Ai( u Oo)(A(u,y) - A(u Q0 ,y))ej for i = 1 , • • • ,ra, j = 1 , • ,n in the sense of 
Corollary 2. Then a is analytic in some neighborhood of 0 and a(0) = 0. 
By Corollary 2 


(4.24) 

I (£ i (A(u) 


A(u 00 ))e i)yvi < 2 

Vv-1 (P + 1} (V + 1} 


2 ( y ; v )^(K 0 T 0 -)w 


From (4.4) 


for l<y + v<m + n. 


(y + v + 1 )(i|> + <|>) 


yv 


(“ + v )a 

v A) ' It 


H 


y+v-2 


(y + 1) (v + 1)‘ 


.2 . , v 2 v v J y+v-1. 


for 2<y+v<m+n, 


since if both y > 0, v > 0 the left-hand side is 0. By Lemma 1 and the 
above 



_< e(y + v + 1) 

Vv-1 



Vv-2 


(V + v )a 

2,.. . ,x2 1 v J y+v-1 


(y + 1) (v + 1) 


2 2 C t V K->y + v> 


(y + 1) (v + l) 


for 2<y + v<m + n 
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Since 


( <+ + *\) vv \ i a ° 1£ “ + 


we obtain 


(<♦ + *\) 


yv 


2e 


H 


(y + l) 2 (v + l) 2 V 


[y + V )(z a + za n ) , 
'» J v 0 y+v’ 


(4.25) 


y+v-1 


for l<y+v<m+n. 


Remark : If y + v = 1, say y = 1 and v = 0, then (4.25) is simply 


N). 


2e 

- 2 2 a 0’ 


which holds by definition of a n . Let b be a majorant to (u„)B e. for 

u i 00 q j 

all i, j, q. By Corollary 2 


U 


B e.) 


< 


1 


2 ( v ) b ^ K 0 T 0 ^y+v’ 


(4.26) 


1 q j WV ’ Vv-1 (u + 1)2(V + 


for l<y+v<m+n. 


Also |(£. B e.)__| < b_. Next 
1 i q j 00 1 — 0 


l < “« > P.v| H ■ <» + » |<"J Vl.»| 


y+v 


H 


y+v 


< T Q (y + 1) 


1 


(y + 2) (v + 1) 


1 , y + v + 1 ^ 

~T2 l v > 


y+v+1 


(4.27) 


1 (y + v + 1)! ,1 


(y + l) 2 (v + l) 2 vl v! 


L a )y+v* 


for 0<y+v<m+n-l, 
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Similarly 


0 


(4.28) 


| ( W,,v| Vv < (y + i) 2 (v + i) 


1 ( y + v + 1) ! 

2 y! v! '•z 'y+v’ 


for 0<y+v<m+n 


Finally 


9 y q 


< e(y + v + 1) I (u . ) 

H I J H 

y+v y+v-1 


(4.29) 


< e(y + v + 1) 


eT, 


[V + 

2, 2 v J y+v 


(y + 1) (v + 1) 


(y + v + 1) ! 


a 


(y + l) 2 (v + l) 2 v! 111 » +v > 


for l£y + v_<m + n, 


and 


K), 


9y q ' j'o.o 


< V 


We now have 


‘ V - <A <“> " «“oo» + <* + «/ h ) J 

1 1 Vv-l 

'• = *■ * V ) j sup | Z e. | ct 

" (U + 1) 2 <V + l) 2 V ( y ' 1 3' > ,+u 

+ *«0 T 0 “Vv + 2ed ‘l + '»o )T<K 0 T 0 “>)„+»}• 


for 1 <y+v<m+n. 
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where we used (4.23a), (4.24), (4.25), (4.26) and (4.19) and summed over 
q = 1 to d. Also note the left-hand side term is 0 when \i = v = 0. 
Hence by (4.20a) applied to (4.27) and the above, summing over j, we obtain 


l* 


.•(first term in (2.12))1 l < - 

i mn | > , \ 2 y . x 2 

H m+n (m + 1} (n + X) 


(m + n + 1) ! 
m! n! 


(4.29a) 


la 2 + I a. a( K 0 T 0 *) + I a-(za + a Q z).b(K 0 T Q a)! +n 


(m + l)^(n + l)^ 


1>! (q«» * zR(a,z) ) 


m+n+1 


where 


Q(a) = C Q a + a*a(K 0 T Q a) 


R(a,z) = a*(a + a Q )»b(K 0 t Q a). 


Similarly , 


(4.29b) 


^•(second term in (2.12)) 


mn 


m+n 


(m +1) (n + 1) 


_ (m + n + 1)! ( . 

2 m! n! ^ ^m+n+1 


It remains to bound the last term in (2.12). From (4.23a) and (4.23b) 


| ( *n - VuvIh^ i („ + u 2 (v ‘ 0 2 1,1 V V ^Vv. 


for l<y+v<m+n, 
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and |(4>^ - 4 >^)qq| _< 2a q. Hence by (4.19) and (4.26) 


l (< *n - V*i B , e jM„ 


P+V-l 


k: 


~ 7m + i) 2 cv + 1, 2 (“: u )' (a+ v« K o T o“>u, 


for 1 <y+v<m+n 


and K(*„-» 5 )* 1 B q e J ) 00 | <2a 0 b 0 . 

Applying (4,20) to (4.29) and the above and summing, we obtain 


A. ‘(last term in (2.12)) 


mn 


H 


m+n 


(4.30) 


(m + l) 2 (n + l) 2 


(m + n + 1) ! 
m! n! 


((a + a Q )(a + 2a Q )b(K 0 T Q a)) 


m+n 


(m 


(m + 


l) 2 (n 


+ 1 )‘ 


• + m ? n! 1)! (z R ( a (z)»z)) 


m+n+1 


with R = Cq(& + a^Kfl. + 2ap)b. The estimate (4.30) holds for m + n = 0 as 
well by the definition of a Q . Now, (4.29a), (4.29b), and (4.30) together 
yield the second estimate in (4.5). 


Remarks : When m = n = 0 the last term in (2.12) is the only nonzero 

term and (4.30) gives a bound for it. 

Also, we have not considered the terms in the end gaps (2.15) separately 
since they have the same form as (2.12). 
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Lemma 8: Let 


u(s,t,y) = l u m (s,y)t i.e., u n = 0 


m>l 


v(s,t,y) = l v (s,y)t 
m>0 


m 


w(s,t,y) = l w (s,y)t . 
m>0 


Consider the following estimates 


(i) 


u 


o 


a 


m| H (m + l) 2 ” 

m-1 


m > 1 


(ii) 


w < 


1 


m| H , (m + l) 2 “ 

m-1 


m > 1 


(iii) 


T„(m + 1) 

v I < 7- C 

m 'H (m + l ) 2 “ 

in 


m > 0 


with a , b , C > 0 , |w-| < T, b n . Then 

m m m— 0L. — 1 0 

H 0 


(A) If (1) holds for 1 < m < k and (ii) holds for 1 < m < k 


( 4 . 31 ) 


<u») k 


T T K 

< ° - ■ 2 ■ 2 [(^(z))(b(z) + b Q )] k . 


H k-1 (k + X) 


1 then 


If (i) holds for 1 < m < k - 1 then 
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(4.31a) 


2 2 

T Ic 

l ^00 /- 2 / _ %x 

( u )k. ^ « (& (z)) . 

** I TT \ ^ K- 


Vl < k + '>' 


(B) If Case 1: (i) holds for 1 < m < k and (iii) holds for 


0<m<k-l or 


Case 2: (1) holds for 1 m _< k - 1 and Vq = 0 = and (iii) 

holds for 1 < m < k - 1 then 


(4.31b) 


(mO k | < (k + i)(a*c) fc . 


\-l < k + 


(C) _If^ Uj = 0 = and (1) holds for 2 ^ m <C k and (111) holds for 
0 < m < k - 2 then 


(4.31c) 


(uv) 


T T K 

< - - - - 1 - ■ W (k + 1) [a(z).c(z)L . 


k| H k _ 1 ' (k + l) 2 


Proof : The results in A and B follow from Lemma 7 (4.19a), (4.20a) 

respectively by considering only one index, say y _< m, m = k, v = n = 0. 
Part (C) is almost immediate: 


| (uv) kL ^ \ | u m[ . | Vi 

1 H k-1 ra= 2 m-1 


, T o T 1 K o (k + x) S 

■"'k-* 1 (k + l) 2 J 2 amck - m 


T n T. K^(k + 1) 

■ 1 ° o (a«e) k , if a 

(k + i r 


= 0 . 
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Corollary 3: Let u = (u.,»**,u ) with u. (s,t,y) = I u. (s,y)t m and 

m>0 m 

u^ satisfy (i) of Lemma 8 for 1 <, ® <, U* Let a(u,y) = £ a a (u-,y)(u - u_) a 
m a 

be an analytic function in the variables (u,y) and suppose 

|a a (uo(s,y),y)| < a £ , (s,y)efi g x 
where UQ(s,y) = u(s,0,y). Then 

| (a(u(s,t ,y) >y) )jj|^ - ^~777 ( “ (K 0 T 0 ttCz) )) U - 

Proof : 



Remarks : If | a| >_ 2 it follows by (4.31a) that 



x a\ , Cr (V 

V JJ £ 2 

\-l ^ + D 



under the weaker hypothesis that (u^) m satisfies (i) for 1 _< m <_ p - 1. 
Therefore, if a^ = 0 for |a| <2, Corollary 3 is valid under this weaker 
hypothesis. 

We now focus our attention on (2.17), the inhomogeneous term in 
rarefactions. We let &(z) be as in the hypothesis of the Main Majorization 


Lemma. 
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We have from (4.4) 



a 




H 


(4.32) 


U + 1 (l 


(y + l) 


2 '■z 


(7 


for 0 < y < k - 2. 


(4.33) |<+ - V C) u| H i| ( * 

v-l 


- X, t) I < . < r (z<l) 

♦ “Ih _~ v v 2 " _1 “ (P + l ) 2 

\x z 


for 2 < y < k, 


and (* - t) Q = (* “ \ t)j = 0. 
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Similarly 


(4.34) 


I C4> ~ X t) J < - TT (za) , for 2 < y < k. 

1 * y 'H , “ (y + l) Z P ” " 

y-1 


Using (4.31c) on (4.32), (4.33), and (4.34) we obtain 


(4.35) !/(('!' - X t) - (<}. - X t))o ± 


C 0 (k + 1) 2 

\ < — r (a) k' 

t j k (k + D 2 

\-i 


From Corollary 3 


|( V (A(u) - A ( u 0 )) e j) u |^ < * 2 “V 


for 1 < y < k - 1 


2 

with a a majorant for £^(A - A(uQ))e^ and = Kq Tq in this case. Also 


u. - u. \ \ < e(y + 1) /u. - u. \ < C n — + ^ a 

3 3 °hl - ( 3 H H - 0 » 2 u 

rl rl - 

y y-l 


for 1 < y < k - 1 , 


Using (4.31b) Case 2 in Lemma 8 applied to the above two estimates and summing 
over j , we obtain 


c n (k + 1) 

(4.36) |(* i (A(u) - A(u q ) )(u - u Q ) s ) k | < j- (a(C Q a)a) k . 

(k + 1) 
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Note that a(0) may be taken to be zero. 

Since 

A(u) - A(u Q ) - A'(u 0 )(u - u„) = l A (u n )(u - u n )“, 

0 0 0 I | >2 a o 0 


if we let a = £ a. , we then have that a is a majorant for 

^ j>2 J q 


£ i (A(u) - A(uq) - A'(uq)(u - Uq) ) e^ . 


Hence, Corollary 3 and the remark following it imply 


(4.37) 


|*i( A " A(u 0 ) " A '* (u " VK) j 


< ^~2 a O (C 0 a) k 

(k + iy 4 U k 


H. 


k-1 


where Cq here is also a bound for 
Since 


u. 




for 1 < p < k - 1 


and since 


C Q (p + 1) 


(4.37a) |(V u - u 0>s)u| <“ 

' 1 ° S y| H (p + l) 2 y 


for 1 < p < k - 1, 


by Case 2 of (B) in Lemma 8 we obtain 
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(♦t " X + + 8 ( ( +t " V " (<{, t ' V ■ V« 


(4.38) 


. °0 (k + U rn 2, 

— o ^ k' 

(k + 1)^ k 


Corollary 3 implies 


(4.38a) |<‘l » q V # | < 7—7 .2 P (C 0 a) )u for » < 8 < k - 1 


Vl (U + »' 


with b a maiorant for Z. B e. 9 Also 

i q 3 


(4.38b) 


I <* " * > »| H - | (t " w i«L 

P-1 p-2 M 


(P + 1) 


2 (z 


for 2 < p < k. 


Since |(+ - «f>) 1 l < 2a Q , we obtain 


(4.38c) I (4» - <!>)„ I < =■ [(za + 2a n z)l 

1 u, h m _ 1 “ (p + l) 2 0 y 


[(za+ 2a Q z)l for 1 < p _< 


and (i|> - 4 > )q = 0. Also, as in (4.29) 


(4.38d) (u ) < 2__(p+l) a for 1 < p < k - 1, 

% J y H p “ (p + l) 2 y - 


Wlth IsjT ( u j( s » y) ) 0 | < V 
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Applying Case 1 of (B) in Lemma 8 to (4.38d) and (4.38c) we get 





< 


c Q (y + l) 

(v + l ) 2 


[(a + a Q )(z a. + 2a Q z)] y , 


for 1 < y < k. 


Applying Case 1 of (B) in Lemma 8 one more time to (4.38a) and the above, and 
summing over q and j , we obtain 

C Q (k + 1) 

(k + l) 2 

Since using Lemma 1 and then (4.4), 


(4.39) | (Jt ± BC4» - +>“y)k.| < 

\-l 


[b(C 0 a)»(a + a Q )(za + 2a Q z)] fc . 


/ 4> + s(4» - 40 \ I £ e(y - 

\ q q / U ' H . 

y-i 


l) 



J o 


, V" - U 1 . , 

1 5 7 Vl 1 ( U + u 


2 (z<l) u 


for 2 < y < k, 


and since 



+ s(4> - 40 \ 

y q/i 


< c o a o 


we obtain 


(4.40) 



+ s(4> - 40 



< — T (za + a n 

(y + l ) 2 0 


z) 


y’ 


for 1 < y < k 
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Also 


I T e(y + 1) 

feOOj <~ 2~^ + a o\ for 0 <y<k-l. 

ds 2 y 'H (v + D U w 

y 


Applying Lemma 8 , part (B), Case 1 to (4.40) and the above, we obtain 


U + s(i|> 


♦\) ( v. 


H 


C Q (y + 1) 

< o- [(<*• + a n )(z a. + a n z)] , 

“ (y + l) y 


y-l 


for 1 < y < k. 


Applying Lemma 8 one more time to (4.38a) and the above, and summing over q 
and j we obtain 


(4.41) 


I + s(* - <(>) y )u s ) k 

\-l 


c 0 (k +D 

1 5 — [(a + a n ) *(za + a n z)»b(C n ct) ]. - 

(k + IT u u U 


Collecting the results from (4.35), (4.36), (4.37), (4.38), (4.39), (4.41), we 
obtain the first estimate in (4.5). 

We consider (2.22) next. From (4.4) 


I"' ' V.| H ,-77777““ 


for 1 < y < k - 1 


| ((vi - u Q0 ) - (v - v 0Q ) 


2T n 

77777““ 


for 1 < y < k - 1 
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Hence Lemma 8 (4.31a) implies 


(4.42) |* ± [<* t - A + )( ( „ - “ 00 ) - (v - v 00 ))]J < ° 2 (g2) k' 


Next, using Corollary 3 we obtain 


* ± U f o (u ’ y) " f 0 (u 00’ y) “ A(u oo )(u - u oo } ) 


(4.43) 


- (f Q (v,y) - f 0 (v 00 ,y) - A(v oq )(v - v Q0 ))] 




(k + 1) 


— ri — — 

^Oq = Z f 0 z with f 0 a majorant for A^fg at Uqq and Vqq, 


Remarks ; Since f 1 = A the left-hand side of (4.43) is quadratic in 


(u - u nn ). 


Using Corollary 3 we also get 


|( V ( V „, - * q W))„| < CT(c 0 


(f(C n a(z))) t 


for 0 < y < k - 1 


where f is a majorant for ‘l = 1, •••,<!. Since 
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(4.43a) |(* y ) M 


Vl (P + x) 


2 (za+ a 0 z) y 


for 1 < y < k, (see 4.40) 


by Lemma 8 (4.31) 


(4.44) |(* 1 (f(u) - f(v))** y ) k | < (f(c Q a(z)).( z a+ a Q z)) k . 

H k _^ (k + 1) 


Putting (4.42), (4.43), and (4.44) together we obtain the third estimate in 
(4.5). 

The expressions for L ^ , and L p from (2.25a) and (2.29) respectively 
are similar so it suffices to restrict our attention to L p in (2.29). Let 

V"-y-V * l . (A p ) a,j (u 00’ y ' 0)(u - »00>“ *y- Then 

01 » J 


(4.45) L 


P <U>y 'V ‘ <V«,0 (u - “00>° + i <X p ) a, j - (u - "00>“ *y' 


j=o 


j>l 


Let X p (z,w) be a majorant of X p , that is 

|„|. t l (x p ) o,j( u oo (y) - y )l l <x P h,j’ y£ V Slnce 


J 


N>p|„ - ... °,,2 t (zat a o z)J Jp for 1 < u < k - 


1 Vl (M + l > 


(see (4.43a) and Lemma 8 (4.31a)) and by Lemma 8 (4.31a) 
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— 2 

(p + i r 



) 


p 


for 0 < y < k - 1 


we obtain by Lemma 8 (4.31) 


|« u " u oo )<X 



c j c' a l 

<5o5> (a H ( ,a + 

(k + i r 


Hence we may estimate the terms in (4*45) by 


(J, " U 00>“ +y)k 


H. 


k-1 


< 1 2 I <Vt i °0 + 

(k + i r j>i p ,J 

£>0 


(k + 1) 


2 zC 0 


j | 0 <V«J« c 0 c > 

£>0 


— o (zR(a,z) ),. 


(k + 1)‘ 


Finally, we estimate the first term in (4.45): 



a >2 


<X p ) a,0 -(u - W 




1 

(k + 1) : 


0 


a| 2 0, 


a o 2 )j )k- 


a o z);l \l 

*(zfl + Sq z)^ 




The fourth and last estimate in (4.5) follows. 
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